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ABSTRACT 

Context. Excess levels of '^N isotopes which have been detected in primitive solar system materials are explained as a remnant of 
interstellar chemistry which took place in regions of the protosolar nebula. 

Aims. Chemical models of nitrogen fractionation in cold clouds predict an enhancement in the gas-phase abundance of '''N-bearing 
molecules, thus we have searched for '^N variants of the NiH^ ion in L1544, which is one of the best candidate sources for detection 
owing to its low central core temperature and high CO depletion. 

Methods. With the IRAM 30 m telescope we have obtained deep integrations of the N'^NH^ (1 - 0) line at 91.2 GHz. 
Results. The N'^NH^ (1-0) line has been detected toward the dust emission peakof L1544. The '''N/'''N abundance ratio in N'''NH^ 
resulted 446 + 71, very close to the protosolar value of ~ 450, higher than the terrestrial ratio of ~ 270, and significantly lower than 
the lower limit in L1544 found by Gerin et al. (2009, ApJ, 570, LlOl) in the same object using ammonia isotopologues. 
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Q-i 1. Introduction 

O ■ The isotopic '"^N/'^N ratio has been measured in a variety of solar system bodies, including giant and rocky plan ets, comets, 
and meteorit es. Its value e xhibits large variations, depending on the selected object and the source for the measurement ( Owen et"al] 
2007h . with the global Jovian value of 450 +100 representing probably the best approximation of the protosolar 



I i ' value (Ipouchet et aL l2004l) . This value differs significantly from the terrestrial one of ~ 272, and even larger ^^N excess was found 



in primitive solar system materials like meteorites, interplanetary dust particles (IDPs), and cometary dust particles returned by the 
Stardust rmssion (Alexander et al. 1998; Messenger 2000; Messenger et al. 2003; Aleon & Robert 2004; Clayton & Nittler 200^ 
> ' iMcKeegan et al.l l2006); to date, the largest '^N enhancement detected is a ''*N/'''N ratio of 65 + 1 1 (corresponding to about four 
. times the terrestrial abundance of '^N) in the hotspots of the meteorite Bells (.Busemann et al.ll2006) . The lack of any significant 
'^C enhancement in the material with the largest ^^N content (iFloss et alJl2"004l) ' rules out a nucleosynthetic origin for the nitrogen 
fractionation and suggests that, as with deuterium, the '^N-rich material results from low-temperature gas-phase ion-molecule 
. reactions and catalysis on cold interstellar dust grains (e.g. Charnlev & Rodgers 2002). 
^""^ The early model of '^N fractionation in the interstellar medium (ISM) developed by Terzieva & Herbst (2000) p redicte d only a 
modest enrichment and was thus unable to account for such large enhancements. Subsequently. , Charnlev & Rodeersl (l2002h showed 
that much higher '^N/'^^N ratios can be generated in high-density cor es, where CO is depleted onto dust grains, but N2 remains 
in the gas ph ase, as appears to be the case in many pre-stellar cores dCasem et al." 1999; Berain et al. 2002; Ta falla etal J 120021: 
iBergin & Tafalla.2007.) . The key fractionation processes are the exothermical reactions (.Terzieva & Herbst.2000.) 
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i^N+'^NzH^ ±^ '^N+'^N'^NH^; (I) 

i^N-t-''*N2H+ ±^ ''*N+''*N'5NH+; (2) 

which, at lo w temperatures, drive '^N into molecular nitrogen through the dissociative recombination of dyazenilium ions 
dMoIek et al.l i2007). Under normal interstellar conditions, '"*N and '^N atoms are continuously exchanged into molecular nitro- 
gen through the sequence 

He+ OH N 

N2 > N — > NO N2 , (3) 

but in heavily depleted regions, there is insufficient OH to drive the above sequence. In these conditions, '^N is preferentially 
incorporated into gas-phase N2 through the dissociative recombination of '^NNH^ and N'^NH^, and into solid NH3 through the 
production of ^^N^ via ^^N N + He^, successive hydrogen ation of ^^N"^, production of '^NHs and freeze-out onto dust grain surfaces 
(ICharnlev & Rodgers 2002: Rodgers & Charnlev 2008b). Assuming selective depletion, this chemistry leads to accretion of ammo- 
nia ice with high '^N enrichment, up to one order of magnitude with respect to the elemental '^N/'"^N ratio, more than sufficient to 
explain the largest measured enhancements. However, this theoretical model is still to some extent speculative, due to the paucity 
of observational data on key '^N-bearing molecules in dense prestellar cores. In this context, the observation of '^N-containing 
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diazenylium ion in selectively depleted dark clouds is an efficient probe to assess if the nitrogen fractionation process is at work in 

the way described, as the model predicts an enhancement in the abundance of both '^NNH^ and N'^ NH'^. 

Previous searches of isotopic variants of NiH^ have been carried out a long time ago (Womack et al.|[T992t iLinke et al.lll983h 
and the detection was successful only toward massive star forming regions, owing to the low sensitivity achieved and also because 
of the selection of sources, which did not concentrate on heavily CO-depleted, centrally concentrated cores (not known at that time). 
To date, no observational data on '''NNH^ and N'^'NH^ in the ISM are available, thus we initiated a survey of '^N-diazenylium in 
cold quiescent clouds, starting from LI 544, which we expected to be a very good candidate source for the detection because (/) the 
deuterium enhancement in this source is very large (N 2D^/N2H^ ~ 0.25, jCrapsi et al. 2007); (//) deuterated species are excellent 
tracers of the high density gas in the center of the core dCaselli et al.ll2002al) . where CO is more heavily depleted and thus where the 
highest '^N enhancement i s expected from th e above reasoning; and (Hi) recent observations revealed that its central region has a 
temperature of only 5.5 K jCrapsi et al.ll2007l) . thus, due to the small zero-point energy changes associated with '^N-fractionation, 
these low temperatures are expected to yield higher '^N/'^N ratios. We report here the positive detection of the N'^NH^ (1-0) 
emission at 91 .2 GHz in this cold dense molecular core. 



2. Observations 

Th e millimetre and submillimetre spectra of N'^NH^ and '^NNH^ were recently investigated in the laboratory bv lDore et al.l 
(l2009l) . and the data were used by the authors to pro duce a ccurate hyperfine line lists adopting the quadrupole coupling and spin- 
rotation constant of the parent sp ecies ((Ca selli et aL Il995h . The same data have also been included in the Cologne Database for 



Molecular Spectroscopy (CDMS, iMtiUer et al.li2001 



l2005h . where lists of hyperfine-free rotational transitions of '^N-containing 



isotopologues of the dyazenylium ion are also presented. 

The observations toward the quiescent Taurus starless core L1544 were performed with the IRAM 30 m antenna, located at 
Pico Veleta (Spain) during one observing session in June 2009. Since the (1-0) transitions of N'^NH^ and '^NNH^ have rest 
frequencies of 91205.6952 MHz and 90263.8360MHz, respectively, it was not possible to observe both lines simultaneously with 
the same detector settings due to the current limitations in the telescope hardware. In the initial observing strategy, a splitting 
of the telescope time into two separate runs was planned; but owing to the unstable weather consitions, we decided to employ 
the whole allocated time integrating the N'^NH^ (1 - 0) transition in order to obtain a spectrum with sufficient signal-to-noise 
ratio. We used the EMIR receivers in the E090 configuration, observing the N'^NH^ (1-0) line in the lower-inner sideband. The 
observations were performed in frequency switching mode, with a throw of +7 MHz; the backend used was the VESPA correlator 
set to a spectral resolution of 20 kHz (corresponding to 0.065 km s ') and spectral bandpass of 20 MHz. Telescope pointing was 
checked every two hours on nearby planets and bright radio quasars and was found accurate to ~ 4"; th e half power beam w idth 
(HPBW) was 27". Scans were taken toward the peak of the 1.3 mm continuum dust emission of L1544 jCaselli et al.ll2002al) . the 
adopted coordinates were RA(2000) = 05''04'"17.2P, Dec(2000) = 25°10'42.8". We integrated for a total of 27.25 hours, with 
two orthogonal polarizations simutaneously observed and averaged together to produce the final spectrum. The rms noise level 
achieved was about 2.5 mK, allowing for a clear detection of N'^NH^ (1 - 0) emission line toward L1544, as illustrated by Fig.[T] 
The spectrum is presented in units of Tmb and was corrected assuming a source filling factor of unity and using the forward and 
main beam efficiencies appropriate for 91 GHz, F^ff = 0.95 and Beff = 0.75, respectively. 



3. Results 

The data processing was done with the GILDASQ software (e.g. I PetvllIOOl : due to the wavy background produced by the 
frequency switching observing method, extensive polynomial baseline subtraction had to be applied to obtain reasonable flat spectra. 
Since the nuclear spin of i^N is i, N'-**NH+ has only one quadrupolar nucleus, '"^N with 7=1, thus its (1 - 0) rotational lines are 
split into a triplet, m aking its detection easier than of the parent species whose hyperfine structure is spread over seven components 
(ICasenietal.lll995h . 

Figure [U displays the averaged spectrum taken toward LI 544. The two stronger F = 2 - I and F = I - I transitions are 
clearly seen, while the weak F - - 1 component is detected at 2cr level. Average line parameters can be estimated by fitting 
Gaussian profiles to the detected lines with the HFS routine implemented in CLASS, which allows to take into account th e hyperfine 
compo nents self-consistently. Adopting the hyperfine splittings and intensities of the 7=1-0 transition calculated bv lDore et al.l 
( 120091) . the HFS fit gives a systemic velocity Vlsr = 7.299 + 0.017 km s"' and an intrinsic line width Av = 0.533 + 0.054km s~'; 
it also indicates, as expected, a low optical depth (t < 0.1) for the line, thus no information on the excitation temperature Tex can 
be derived from the observational data. The results of the HFS fit are summarised in Table [1] and the resulting spectral profile is 
superimposed in Fig.[T]as a dotted trace. 

The column density of N'^NH^ has been calculated from the integrated line intensity of the strongest F = 2 - I component, 
which exhibits the best spectral profile. The Gaussian fit gives f Tnibdv = 14.2+ 1.1 mKkms-i and Av = 0.575 + 0.055kms-i. 
We adopted the constant excitation temperature approximation (ICaselli et al. 2002b), assuming a Tex val ue of 5 .OK as derived from 
observations of the N2H^ (1-0) hyperfine structure toward the same object and the same offset position (.Caselh et al..,2002a) . From 
the solution of the radiat ive transfer equation with the assumption of optically thin emission, one has the following expression for 
the total column density dCaselh et alJl2002bl) 

ex l±-]-lY\^ /_E^\^/Xbdv_ 

.^""^ I fc^ex / J ^''P 1^ kT,, I 7,(rex) - Jy(T^g) ' 



87rv^ 6sr(7^ex) 

c^A 2Fu + 1 



' See GILDAS home page at the URL: http : //www . ir am . f r/IRAMFR/GILDAS. 
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where A is the emission Einstein's coefficient for the hyperfine transition, Jv(T) is the radiation temperature of a black body at 
temperature T, is the upper state energy, and 2,si(rex) is the spin-rotational partition function for the excitation temperature 
obtained by summing over all energy levels of importance. 

The A coefficient for the F — 2 - \ hyperfine line was calculated from the formula given bv iDore e t al. (2009) and with the 
weighted average of the literature values of the N'^'NH^ dipole moments, yu - 3.31 + 0.20D dHavenith et al..,1990 . Table II), 
resulting A - 3.23 + 0.40 • 10"^ s"'. Substitution of all terms of Eq. (|4]i gives a total column density of N'^NH^ toward L1544 of 
(4.1 + 0.5) ■ 10'" cm"^, where the estimated uncertainty is obtained by propagating the errors on A, Tgx, and the integrated intensity 
derived from the present observations. 

4. Discussion 

The column density of the main dyazeniulim ion derived toward the "dust peak" of LI 544 bv ICrapsi et al.l (l2005h is (18.3 + 
1.9) ■ 10'^ cm"^ and thus the resulting [NaH^/N'^NH^] abundance ratio is 446 + 71, which is well compar able with the recognised 
protosolar value of the '"^N/'''N ratio, as measured in th e Jupiter atmosphere (~ 450, iFouchet et al.ll2004l) . or in osbornite-bearing 
inclusions from meteor ites (~ 420, Meibo m et al.ll2007l) . 

Very recently Gerin et alj (^09) reported on a search for '^NHaD in dense cores, with the aim of measuring the nitrogen isotopic 
ratio in the ISM. They observed several sources, including L1544, obtaining '"^N/'^N ratios ranging from 350 and 850. In L1544 
the detection of '^NH2D was not achieved, leading to an estimation of a lower limit of 700 for this ratio. The value we obtained for 
N2H^ is nearly two times smaller than this estimate, indicating that the mechanism of nitrogen fractionation at work in these cold 
dense cores produces marked differences of '^N enhancement among differen t chemical species. 

The time dependent coupled gas/solid chemical model of ICharnlev & Rod gers (2002) (see also iRodgers & Charnl"evll2008bl) 
predicts that at the end of evolution significant amounts of '^N-rich ammonia are frozen onto ice mantles, while the gas phase 
becomes enriche d at early times, before the complete depletion of molecules. Improved models dRodgers & CharnlevI l2008al : 
iGerin et ak I l2009h which include additional ion/neutral and neutral/neutral reaction channels predict that assuming typical dense- 
core para meters, '^N enrich ment of ammonia is only moderate in the gas phase, while much stronger enrichment is expected 
for N2H^ (iGerin et al ] |2009[ figure 2). Thus, it appears that our finding is not compatible with th e above picture, given that the 
[N2H"^]/[N'^NH'^] ratio in L1544 is consistent with the '''N/^^N abundance ratio in the local ISM d Wilson & Roodlll994 . In any 
case, further observations are needed to effectively test these chemical models. In particular, the detection of the other isotopologue 
'^NNH^ will be very interesting, since due to the different exothermicities of the reactions ([T]i and (O, the two ions should be 
fractionated to a different degree. IRodgers & CharnlevI (i2004i) predicted a gas-phase ratio [N'^NH^]/['^NNH^] of about 2, which 
differs significantly fro m the equal abund ances that would exist without the fractionation mechanism. A tentative determination 
of this ratio is found in lLinke et al.l (Il983h who obtained 1.25 in DR21 (OH), but it is likely that new determinations in centrally 
CO-depeleted cores might yield higher values. 
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Table 1. Results of the hyperfine fit on the observed spectral profile of the N'^NH^ (1 
parentheses refer to Icr uncertainties in units of the last quoted digit. 



0) transition toward L1544. Numbers in 



Line 

F' -F 


Rest frequency" 
(MHz) 


A coefficient 
10^ 


Vlsr 
(km s-') 


(mKkm s"') 


Av 
(km s-') 


lOV* 


0-1 


91208.5162 


3.23(40) 


-1.97F 


2.68(34) 




2.31(20) 


2-1 


91205.9908 


3.23(40) 


6.329(17) 


13.4(17) 


0.533(54)'' 


11.6(11) 


1 - 1 


91 204.2602 


3.23(40) 


12.018" 


8.1(10) 




6.95(63) 



" From Dore et al. (2009). 
* Optical depth derived assuming r„ 



5.0: 



0. 1 K from lCaselli eTH] ( l2002ij) . 

" The hyperfine splitting with respect to F = 2 - 1 component was kept fixed in the HFS analysis. 
Gaussian FWHM. Assumed equal for all the components. 
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Fig. 1. Spectrum of the N'^NH^(1 - 0) transition observed toward L1544 (grey trace), and computed spectral profile resulting 
form the HFS fits (red dashed trace). The superimposed histogram (black trace) indicates the position and relative intensity of the 
hyperfine components. 



